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Critical Behaviour near the Mott Transition 
in a Two-Band Hubbard Model 
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The Mott metal-insulator transition in the two-band Hubbard model in infinite dimensions is 
studied by using the linearized dynamical mean-field theory. The discontinuity in the chemical 
potential for the change from hole to electron doping is calculated analytically as a function 
of the on-site Coulomb interaction U, and the charge-transfer energy A between the d- and 
p-orbitals, transfer integrals t p d, t pp , tdd between p-d, p-p and d-d sites respectively. The critical 
behaviour of the quasiparticle weight is also obtained analytically. 
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§1. Introduction 

The Mott metal— insulator transition (MIT) is ob- 
served in various transition metal oxides and have been 
extensively studied for years. The dynamical mean-field 
theory (DMFT), 1,2 ) which becomes exact in the infinite 
spatial dimensions, has been used in various models and 
have given plausible results for the MIT. However, it is 
difficult to obtain the critical behaviour of the MIT nu- 
merically as the quasiparticle peak becomes extremely 
sharp in the limit of the MIT. There is an alternative 
way to resolve the adifficulty. The linearized DMFT 
(LDMF) 3 ~ 6 ) is an approximate but analytical method 
which becomes the most correct on the MIT point. 

There are two types of the MIT in real materials: 
the Mott-Hubbard (MH) type in which the energy gap 
is given roughly by the Coulomb interaction and the 
charge-transfer (CT) type in which the energy gap is 
given roughly by the charge transfer energy. In order to 
describe the two types of the MIT, we need at least the 
two-band Hubbard model which is characterized by the 
Coulomb interaction U, charge transfer energy A and 
p-d, p-p and d-d transfer integrals t pd , t pp and tdd- 5 ^ 

The critical behaviour of the MIT has been studied 
using the LDMF for the single-band Hubbard model. 3, 6 ' 
It has also been studied for the two-band Hubbard model 
in our previous paper. 4 ) However, effects of t pp and tdd 
were not considered there, which are not negligible in 
actual compounds but make important contribution to 
the critical behaviour near the MIT, as found to be sig- 
nificant to determine the phase boundary of the MIT. 5 ) 
In the present paper, we study the discontinuity in the 
chemical potential and the critical behaviour of the MIT 
in the two-band Hubbard model in the presence of t pp 
and tdd- 
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§2. Formulation 

We consider the two-band Hubbard model, 5 ) 

vq (id),* 

^ (i,i'),<7 i,cr 

+ yE (pipy- + h-c) + «*>$>W' ( 2J ) 

where d\ a and pj a are creation operators for an electron 
with spin a in the <i-orbital at site i and in the p-orbital 
at site j, respectively. t pd , tdd and t pp are the transfer 
integrals between the nearest neighbor p-d, p-p and d-d 
orbitals, respectively. In eq. (2.1), we assume that p- 
and d-orbitals are on different sub-lattices of a bipartite 
lattice, more explicitly, a Bethe lattice with connectivity 
q. In the limit q — oo, the self-consistency equations for 
the local Green's functions are given by 2 - 1 

^(z)- 1 =z-e d - t 2 pd G p (z) - t 2 dd G d (z), (2.2) 

Gp(z)- 1 = z-e p - t 2 pd G d (z) - t 2 pp G p (z), (2.3) 

where G p (z) is the local Green's function for the p- 
electron and Gd(z) is that for the d-electron; ed = e<jo — 
/x = —fi and e p = e p o — [i = A — jjL, where we set edo = 
and the CT energy A is defined by A = e p — ed > 0. 

In the LDMF, the two-band Hubbard model eq. (2.1) 
is mapped onto a two-site Anderson model 

a 

+ e c ^ctc CT + V £(/jc CT + 4/ ff ), (2.4) 

with e c = and e/ = e f j = — /i. In the model eq. (2.4), the 

non-interacting impurity Green's function is: Goiz^ 1 = 

v 2 

z — ef —.In the limit V — » 0, the local Green's func- 
tions are given by G d (z) = ^ and G p (z) = ^-, 4,5 - ) near 
the Fermi level with small weights Z d — > and Z p — > 0. 
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Then the self-consistency equations (2.2) and (2.3) are 
reduced to a simple equation 



^pd^V + tdd^d — V 2 . 



(2.5) 



To second order in V, the quasiparticle weights for the 
d-electron and the p-electron are given by 7 ' 

Z d = F(U,^)V 2 , (2.6) 

Zp = A(t pd ,tp P ,U,A,n)V 2 , (2.7) 



where 



F(U,f,) = ^+ n ,J + „ (2.. 



A(t pd , t pp , U, A, fj.) 
and 

E p = A - 



2fi 2 2(U-li) 2 Li{U - n)' 



E 2 -t 2 ' 



1 



+ 



1 



(2.9) 



(2.10) 



A - /i pa { 2n 2(U - n) 

Substituting eqs.(2.6) and (2.7) with eqs.(2.8) and (2.9) 
into eq. (2.5) we obtain an equation to determine the 
MIT point, 



P(t p d, t PPl t d d, U, A, jtt) = 1, 



(2.11) 



where 



P = t 2 pd A(t pd , tpp, U, A, M ) + t 2 dd F(U, n). (2.12) 

After self-consistency cycle, V increase exponentially 
with iteration number for P > 1, and then the single 
pole approximation breaks down resulting in the metal- 
lic solution, on the other hand P < 1 correspond to the 
insulating solution. In eq. (2.11), P includes the chemi- 
cal potential fj, which has to be determined explicitly to 
obtain the critical values of the MIT. We can use a cer- 
tain condition to determine [i, based on the fact that at 
the MIT point P has a minimum value as a function of 
/i. 4 ) This condition gives 

g^Pitpd, t PP , tdd, U, A, fx) = 0. (2.13) 

§3. Discontinuity in the Chemical Potential 

By solving eqs. (2.11) and (2.13), the MIT phase 
boundary is determined and shown in Fig. 1 (see Zi/x = 
0). In the insulating regime, the chemical potential fi± 
in the limit n — > 1 ± 0+ is obtained from eq. (2.11): 
P\p=p ± — l, 4 ) which yields the discontinuity in the 
chemical potential, Afi = \i+ — In Fig. 1, we show 
the contour map for A[i on the A-U plane, where the di- 
agonal dotted line represents a boundary separating the 
CT and the MH regimes, and the dashed lines are the 
critical values of U c and A c in the limit of A — * oo and 
U — ► oo, respectively. The discontinuity in the chemical 
potential smoothly connects the Mott- Hubbard type and 
the charge-transfer type insulators. In the limit A — + oo, 
the critical value of U c — 6 t dd is the same as that for 
the corresponding single-band Hubbard model. 3, 6 ) 



§4. Critical Behaviour near the Mott MIT 

We discuss the critical behaviour in the vicinity of the 
Mott MIT at half-filling in the two-band Hubbard model. 
In this case, we need the result of the quasiparticle weight 
for the d-electron up to fourth order in V , 



Z d = F(U,fi)V 2 -G(U,n)V\ 

where F is given in eq. (2.8) and 

29 24 22 

G{U,n) = 



(4.1) 



2^ 4 fi 3 (U-fi) fi 2 (U-fi) 2 
24 29 



(4.2) 



The quasiparticle weight for the p-electron Z p is also cal- 
culated up to fourth order in V, 

Z P =A(t pd , t pp , U, A, fi)V 2 -B(t pd , tpp, U, A, fi)V 4 , (4.3) 

where A is given in eq. (2.9) and 

G .a IF B\Ep 



B(t pd , tpp, U, A,fi) — t 



P d p2 _ f2 
p pp 



+ t: 



pd 



F 2 {2t 2 pp + 3B 2 E p (A-ri} 
+ t p d' 



with 



B x = 



{E 2 p -t 2 ppY{A-ri 
-3 U 3 + 7 U 2 n - 3 U [i 2 + 2 ^ 



tx*(u-nY 



B 2 = l + tt 



1 



1 



(4.4) 

(4.5) 
(4.6) 



pd \2fi 2 ' 2{U-n¥ J ' (A-m) 5 

Substituting V 2 from the self-consistency equation (2.5) 
into eqs. (4.1) and (4.3), we obtain 

P- 1 



Z d = F(U,iJ)- 



Q 



(4.7) 
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Fig. 1. The contour map for the discontinuity in the chemical po- 



tential An = 0, 1, 2, 3, 4, 5, 6, 7 at t 



pd — tpp 



1 and t dd = 0.2. 



The dotted line represents a boundary separating the MH and 
CT regime. 
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where 



In the limit U c — > oo, the coefficients are 



Q = t 2 pd B(t pd , t„, U, A, fi) + t 2 dd G(U, fi). (4.8) 

At the MIT point with the critical values U c , A c and 
fi c , eqs. (4.7) yield Z d = from eq. (2.11). When A or U 
decreases from the MIT point, the quasiparticle weight 
for the d-electron Z d eq. (4.7) increases as 



C A = 1.31, C v = 0, 



(4.11) 



Z r i — Ca 



A 

U 

1 



(4.9) 



(4.10) 



near the MIT point at half-filling, respectively. The co- 
efficients Ca and Cjj are shown in Figs. 2(a) and (b). 

When t pp increases, the p-component of the quasiparti- 
cle weight increases, then, Ca increases and Cjj decreases 
(compare the dotted line with the solid line in Figs. 2(a) 
and (b)). On the other hand, when tdd increases, the 
d-component of the quasiparticle weight increases, then, 
Ca decreases and C\j increases (compare the dashed line 
with the solid line in Figs. 2(a) and (b)). 

In the CT regime (U > A), the MIT occurs at A = A c 
when A is varied for a fixed U. As A decreases below 
A c for a fixed U, the quasiparticle weight increases as 
given in eq. (4.9). With increasing U c (decreasing A c ), 
the coefficients both Ca and Cu decrease due to the 
increasing correlation effect as seen in Figs. 2(a) and (b). 




Fig. 2. The coefficients Ca (thin lines) and Cjj (thick lines) of 
the quasiparticle weights for the ci-electrons Z d = Ca (l — ■£-) 
and Z d = Cjj (l — jj-^j , respectively, as functions of A c (a) and 
U c (b) for several values of t pp and t d d with t pd = 1. The thin 
vertical dotted line shows a boundary of MH and CT for t pp = 1 
and t dd = 0.2. 



for t pd = tpp = 1 and t d d = 0.2. 

In the MH regime [U < A), the MIT occurs at U = U c 
when U is varied for a fixed A. As U decreases below U c 
for a fixed A, the quasiparticle weight increases as given 
in eq. (4.10). With increasing A c (decreasing U c ), Cjj 
increases due to the decreasing correlation effect, while 
Ca decrease due to the decreasing p-component of the 
quasiparticle weight (See Figs. 2(a) and (b)). In the limit 
A c — > oo, the coefficients are 



C A = 0, 



Cr 



(4.12) 
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for tpd — tpp — 1 and tdd — 0.2. In the limit A c — > oo, 
the effect of the p-band becomes irrelevant and, then, 
Ca = and C\j = jf is the same value as that for the 
single-band Hubbard model. 3 ) We note that, even in 
this limit, Ca is finite and Cjj ^ j| in the case with 
tdd = 0, because the electron has to transfer between d 
and p orbitals through the hopping integral t p d and the 
effect of the p-band is still relevant as shown in Ref. 4. 

§5. Discussion 

In the real materials, the ratio of the transfer inte- 
grals are estimated as: La2Cu04, t p d ■ t pp = 1 : 0.50, 8 ) 
LaTi0 3 , tpd : tpp = I : 0.15, LaNi0 3 , t pd : t pp = 1 : 
0.27 9 ' 10) (no data for tdd)- There is little data for the 
critical behaviour of the MIT when U or A is varied 
at the half-filling because of the difficulty of the exper- 
iments. n ) Further progress in both theoretical and ex- 
perimental study for the critical behaviour of the MIT is 
expected. 
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